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Abstract 

It will be shown how to map a simple one-dimensional tight binding model with a cosine po- 
tential in one dimension exactly to a two dimensional tight binding model with periodic boundary 
conditions with the presence of a single flux quantum spread evenly on the torus. The mapping 
is is achieved by a partial sequence of "Fast Fourier Transform" (FFT) steps which if completed 
would be an exact Fourier transform of the original model. Each step of the FFT recursively maps 
a tight binding model into two decoupled sublattices of half the lattice length. 
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I. MOTIVATION 



The following paper is dedicated to Nihat Berker and describes a mathematical curiosity 
that I stumbled across a few years ago. It brings together two ideas that were a daily topic 
of discussion with Nihat during my graduate student days: duality and renormalization. 
The work also resulted from my admiration of the elegant recursive algorithms that lie at 
the basis of the efficient solution of a number of practical problems, most notably sorting 
and the Fourier transform. 

For a transitionally invariant problem, a Fourier transform usually either solves or at 
least decouples the spatial degrees of freedom so that the original problem becomes simple. 
When applied directly on data of length it uses a number of operations that scale like 
A^^. The "Fast Fourier Transform" (FFT) dramatically reduces the number of arithmetic 
operations using a clever recursive algorithm. An interesting question would be if if a single 
FFT recursion could somehow be the basis of a renormalization transformation that would 
partially decouple a problem at each stage of recursion. Ultimately I was not able to answer 
this question but in the process came up with a surprising result; that repeated application 
of the FFT recursion can map a problem from one to two dimensions. 

II. BACKGROUND 

Before proceeding to the calculation this I make a brief comment on the history of the 
FFT. In 1965 Cooley and Tukey wrote a famous paper entitled "An Algorithm for the 
Machine Calculation of Complex Fourier Series" . pLj Cooley and Tukey cited I.J. Good [2] 
as having influenced their work. The impact of the "Cooley- Tukey algorithm" was enormous 
and their "Fast Fourier Transform" immediately became indispensable tools for engineers 
and scientists as computers started to become widely available at this time. These days 
the "FFT" algorithm is part of every numerical package and most of us don't pay much 
attention to the algorithm. 

Soon after the paper came out it was realized that the "FFT" algorithm had been re- 
peatedly rediscovered earlier. A detailed historical discussion on the subject is available by 
Heideman, Johnson and Burrus.[3] The "FFT" was eventually traced back to C.F. Gauss. It 
appears certain that he discovered the algorithm around 1807 but did not publish it. Given 
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the impossibility of making significant numerical calculations at the time Gauss probably 
did not consider the result particularly important. The algorithm was eventually published 
posthumously in his collected works. [1] Thus Gauss' original work was done approximately 
simultaneously with Fourier's development of his famous trigonometric expansion. In my 
opinion though, Cooley and Tukey deserve great credit for realizing in 1965 the importance 
of an efficient numerical Fourier transform and deriving the algorithm at a time when it 
became useful. 

The power of the FFT algorithm is that it dramatically reduces the number of arithmetic 
operations required for a Fourier transform of a large data set. For the most common "radix- 
2" algorithm the number of data points is taken to be L = 2". The number of operations will 
be on the order of 2nL rather than the number of operations that results from applying a 
simple matrix operation. The FFT algorithm can also be implemented on arbitrary length 
data where the number of arithmetic operations scales like L times the sum of prime factors 
of L. When L is prime, the FFT reduces to the ordinary matrix operation. 

Without going into detail, I will now write down the key formulas using radix-2 algorithm 
that shows how the FFT is based on recursion. We begin with data Xm where < m < L 
and L = 2^. We define the Fourier transform 

V -tv m=0 

The sum can be split into two parts: 

^ L/2-1 L/2-1 

Vk = -/kK 1^ a;2me^/2 +6^ ^ 0:2^+16^2 ) 

v2 rn=0 m=0 

27ri 7, 

= Ck + e L Ok 

where can be seen to be the Fourier transform of the even indexed subset of and is 
the Fourier transform of the odd subset. This procedure can therefore be done recursively, 
splitting Cfc and Ok successively into into odd and even subsets each with half as many data 
points. 

This notation is clumsy and does not lend itself to recursion. To implement it systemati- 
cally we must keep track of which order of even and odd blocks we use to construct Xk as we 
apply the algorithm. We use the label Cj^k keep track of our recursive indexing. We begin by 
defining cj^i = xj to consist of the original data. The recursion then consists of breaking the 
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first index j into odd and even blocks. Since I will shortly write down a precise formula, I 
will first just the sketch the idea of constructing the indexing. We lay the "odd" and "even" 
blocks next to one another and put in a binary index "0" last if the block is even and "1" 
if the block is odd. The original block is periodic in length L; each new block is periodic in 
length L/2 and the number of blocks, say in number is doubled to length 2N. Thus the 
original index {j, 1) will be replaced by an index {j', k) where the "block" index k encodes 
odd and even block and (j') describes the new order in the subblock. Doing this recursively, 
the index k will gradually become the "frequency" or "momentum" index. The recursion 
is finished when the entire space index has been reduced to length one and the momentum 
index has become length L. [5] 

Let us use the notation r : L to denote the integer r modulo L. We encode the algorithm 
as follows: 

Cr:L,s:N = {cr:L/2 , s:2N + G^""" ^ ^ Cr:L /2 ,{s+N):2n) (l) 

The inverse is given by 

1 J {.Cr:2L,s:N/2 + C{^r'+L):2L , s:N /2) if S<N/2 

Cr:L,s:N — < _ . , {^) 

V2 '''''1^ {Cr:2L,s-N/2:N/2 " C(^r+L):2L , s-N/2:N/2) if S>N/2 

It is straightforward to apply this operation to quantum mechanical operators Cr:L,s:N 
where c is simply the one particle (destruction) operator on an L x lattice of "pseudo 
real space" dimension L and "pseudo k-space" dimension A^. For the initial iteration L is 
the original space dimension, and = 1. Equivalently we can simply interpret as the 
column vector being transformed and cj as a row vector. 

Thus we begin with a periodic "real space" lattice of length 2" to which we initially attach 
a bogus "momentum space" lattice of length 1 to begin the recursion with Cj^i. At each FFT 
step we halve the length of the real space and double the length of the momentum space. 
Eventually when we are done applying the recursion n times the "real space" lattice will 
have shrunk to length one and the variables Ci^k thai correspond to the momentum space 
lattice will be the Fourier transform of the original variables. 

What if we stop this recursion "half way" ? Provided L is 2^*" where m is integer, we will 
obtain a two dimensional periodic 2™ x 2™ lattice index when we stop the procedure after 
m recursions and thus will have obtained a mapping to a "two-dimensional" 2"^ x 2"^ index 
Cj^k on a torus from our original 2^™ periodic index. 
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Clearly this transformation is highly nonlocal since it mixes Fourier components of indexes 
very far apart. If we start with an arbitrary real space operator written in terms of the real 
space operators and apply this prescription it is tempting to conclude that we will always 
get a terrible mess during the recursion. This is indeed true generically. However, we 
know that a transitionally invariant operator is diagonalized by the Fourier transform, i.e. 
written as a simple local operator in the fully transformed basis. Thus eventually after 
recursion the momentum space blocks must decouple and the intuition that the problem 
gets more and more intractable under recursion is violated. Furthermore, a model which is 
invariant under Fourier transform should be symmetric in momentum and space coordinates 
half way through the recursion, so perhaps such a model will retain a structure during the 
intermediate steps of the recursion as well. It is precisely question we will explore further 
by example. 

In the rest of this communication, I will focus on the following model under the "partial" 
FFT operation. We begin with 

N-l N~l 
hNxl = (4cr+l:Ar + Cc) - 2 ^ 4c,. COS 2'Kr/N. (3) 

r=0 r=0 

This is in fact a special case of "Harper's model" where periodicity is exactly one relative 
to the length of the lattice, p] a point which I will discuss at the end of this paper. 

We know that the first term transforms to the cosine under Fourier transform and vice 
versa. For this choice of parameters the model is "self dual" i.e. transforms to itself under 
Fourier transform. It can therefore be expected to be symmetric in some way when partially 
transformed n times, where 2^" = L. For illustration, we begin with a 64 x 1 lattice and 
apply the FFT step three times to transform this model numerically to an 8 x 8 lattice. The 
result of applying this transformation is given by the following Hamiltonian 

/i8x8 = - XI U,Sclcr+5. (4) 

r,5 

where 5 are only nearest neighbor lattice vectors, where \tr,5\ = t and the phases of t 
correspond to a single flux evenly spread throughout the entire lattice. [7] The value of 
try = e'^™r,r'/^'^ where Ury is shown in Fig. [ij 

The vertical bonds with phase zero represent the mapping of the cosine potential, whereas 
the horizontal bonds that carry the flux is the transformation of the hopping. 
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Let us look more carefully at the distribution of fluxes. With the single exception of a 
plaquette in the top row, we note that the sum of integers around each plaquette is one, so 
that the product of phases around each plaquette is exactly e^'^'^l^^ which I interpret as a 
single flux quantum evenly spread on the torus. Periodic boundary conditions are fulfilled 
and the left and right boundary are mapped to each other, likewise for top and bottom. The 
rule that the sum of phases around each individual plaquette is the sum of phases around the 
entire boundary is obtained form the fact that a plaquette on the top row has flux —63 ~ 1. 
An equivalent phase choice of 2-7? instead of zero for the top link is equivalent to assigning a 
flux of 27r/64 through each plaquette. 

Although I do not present further results here, I assert that similarly applying the FFT 
recursion a different number of times result in rectangular lattices all with a total 27r flux 
threading each periodic lattice. 

It is interesting to compare the eigenstates of the tight binding model to the corresponding 
states mapped by the partial FFT. The absolute square of the six lowest eigenstates ■i/'j of 
of /i256xi is plotted in Fig. |2} These corresponding eigenstates of /iig is plotted in Fig. [3] 
and thus represents the "2-D" mapping of the "1-D" states in Fig. |2j The calculations are 
displayed here for values of which are not large; analogous results for larger values of A^ 
have been computed. 

Connection to Harper's equation 

It is well known that an infinite two-dimensional tight binding model in a constant mag- 
netic field maps to a one tight binding model with a periodic or quasiperiodic potential. In 
addition to P.G. Harper's paper mentioned previously, there are several other early classic 
papers on listed in the references.^ 

In 1976, R. Hofstadter's performed a numerical analysis of the problem and pointed out 
the beautiful self-similar structure of the energy spectrum. [9] A number of paper followed 
in the 80 's that explored this using renormalization group ideas. [10] The extent to which 
these are pertinent to the present problem is addressed in this section. 

For a tight-binding model in the presence of a magnetic field with flux ajl-n per plaquette 



6 



we can introduce a Landau gauge and write the Hamiltonian as 

h = - tr,5clcr+S (5) 
r,S 

where tr^±y = t and tr,^x = texp{±iTTay). Noting the invariance of the problem for 
translations along the x-direction we introduce plane waves along the x axis and define 
Cmx+ny = e*^™c„ rcsultiug in a Hamiltonian for c„ of the form in 

hiD = — ^(cjcr+i + cc) — 2 ^ clcr COS {2Trra — u). (6) 

r r 

If we demand that this Hamiltonian be periodic with period q we find that that a = p/q 
is the ratio of integers. The flux per square in the 2-D problem Eq. [5] is thus a integral 
multiple of 27r/g. Is this the transformation we have been doing in the previous section? 

If we begin with a p x p periodic two-dimensional lattice and add a constant magnetic 
field described in Landau gauge, we find that u must obey u = 271171/ p with m integer and we 
obtain an equivalent set of p one dimensional problems Eq. [6] with periodicity p indexed by 
the phase u. In a conventional formulation of the Harper problem we therefore cannot reduce 
the flux quantum to sum to a single quantum over the entire 2 — D lattice. Indeed, the flux 
configuration in the 2-D Landau gauge is quite different from that shown in Fig. [3] which 
describes a single flux quantum spread over the lattice and lacks translational invariance in 
both X and y direction. (Note the top row vertical bonds is not uniform along x.) A careful 
analysis of the gauge choice in Fig. [3] or another clever choice of gauge other than Landau 
might possibly be used to show that the FFT remapping of the one-dimensional problem of 
length p"^ is in fact Harper's problem, but I have not been able to do this. 

More digressions on the Harper problem would stray further from the main point of the 
communication which is the presentation of a decoupling scheme of the tight binding model 
that is in principle more general than Harper's problem. Unfortunately I can find no other 
notrivial example than the relation between the one and two dimensional Harper model, but 
with flux and periodicity reduced from the ordinary formulation of that problem. 

Further work 

These proceedings has provided an opportunity to publish work which I find interesting 
but admittedly incomplete. When performing the recursion on the tight binding model 
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only, each FFT recursion decouples the spatial degrees of freedom into two coupled blocks. 
My original hope was that these transformation might be the basis of a renormalization 
transformation; if an interaction could be found that was similarly decoupled, a novel way 
to solve an interacting model would be found. Unfortunately this goal has not been fulfilled 
but perhaps someone might be motivated to look at this problem again with the observations 
presented here. 
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FIG. 1: Integers n^r' where t^r' = e^'^*"''''''^^^ are shown. 
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FIG. 3: Plot of l^/'j./cp for the six lowest energy eigenstates of /iiexie rnapped to the two dimensional 
lattice 
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